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I.  INTRODUCTION 

The  work  described  in  this  report  was  conducted  from  1 October  1974 
to  31  March  1975  under  Advanced  Research  Projects  Agency  Contract  DAAB  07- 
74-C-0470,  P.M.  Raccah,  principal  investigator,  F.  Rothwarf,  and  G.  Iafrate, 
contract  monitors.  The  purpose  of  this  investigation  is  to  elucidate  the 
physical  mechanism  responsible  for  the  pyroelectric  effect  in  ferro- 
electric crystals.  We  have  proposed  a generalized  molecular  field  theory 

l 2 

(GMFT)  of  ferroelectricity  and  pyroelectricity  * which  yields  good 
agreement  with  the  available  experimental  data  for  BaTiO^.  Unfortunately, 
these  results  are  limited  in  two  respects.  First,  only  BaTiO^  has  been 
studied  in  sufficient  detail  to  provide  the  input  data  needed  for  our 

t 

theory,  and  second,  the  calculations  reported  previously  are  approximate 
ones,  valid  only  when  the  polarization  is  near  zero  or  near  its  saturation 
value.  We  present  here  exact  calculations  based  on  the  generalized  mole- 
cular field  theory"^  which  remove  the  restrictions  on  the  magnitude  of  the 
polarization.  These  were  undertaken  both  to  explore  the  detailed  behavior 
of  the  CM  Pi  and  to  provide  a framework  for  the  analysis  of  measurements 
currently  in  progress  on  a number  of  materials  in  order  to  obtain  input 
data  for  the  theory.  The  experiments  will  be  described  in  a forthcoming 
report . 

After  a brief  summary  of  the  GMFT  and  the  introduction  of  reduced 
(dimensionless)  variables  which  greatly  simplify  the  calculations  (Section  II), 
we  present  in  Section  III  exact  calculations  of  the  static  properties  of 
our  model  (spontaneous  polarization,  critical  temperature,  critical  pola- 
rization). Then  Section  IV  gives  exact  results  for  the  dynamic  properties: 
susceptibility,  pyroelectric  coefficient,  and  response  figure  of  merit 
for  pyroelectric  detection.  Throughout,  we  treat  both  the  case  of  a local 
field  expressed  as  a power  series  in  the  polarization  and  that  of  the  field 
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represented  by  a Pade  approximant^.  We  also  compare  the  exact  results 

with  those  obtained  from  an  approximate  free  energy  of  the  type  proposed 

3 

by  Devonshire  . 

An  important  result  of  the  present  study  (Section  IV)  is  that  the 
GMFT  predicts  that  the  response  figure  of  merit  of  a pyroelectric  detector 
M^,  which  is  the  ratio  of  the  pyroelectric  coefficient  and  the  dielectric 
constant,  will  exhibit  an  interesting,  and  potentially  useful,  universal 
behavior.  When  expressed  in  terms  of  reduced  variables,  this  figure  of 
merit  depends  only  on  the  spontaneous  polarization,  indicating  that 
should  be  regarded  as  a function  of  polarization  rather  than  temperature 
when  trends  or  systematic  variations  between  materials  are  being  sought. 
Work  is  now  in  progress  to  exploit  this  fact  as  a tool  for  optimizing  M . 

Finally,  Section  V presents  a preliminary'  fit  of  the  observed  tempe- 
rature dependence  of  the  spontaneous  polarization  of  triglycine  oulfate  to 
the  prediction  of  the  GMFT.  Good  agreement  is  obtained  by  trial-and-error 
adjustment  of  GMFT  parameters,  and  a computer  program  is  now  being  written 
to  obtain  the  best  fit  (in  the  sense  of  mean-square  deviation)  and  so 
determine  the  GMFT  parameters  with  the  greatest  precision  compatible  with 
the  accuracy  of  the  experimental  data. 
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II.  SUMMARY  OF  GMFT 

As  is  well  known,  molecular  field  theory  predicts  that  a system  of 

■“4 

N independent  electric  dipoles  of  moment  p which  can  assume  any  orientation 
in  space  will  exhibit  a macroscopic  polarization  P at  temperature  T given 
by 

P = Np  (coth  x - l/x)  , (l) 

where 

x = pE^/kT  , (2) 

is  the  local  electric  field  at  the  site  of  each  dipole,  and  k is 
Boltzmann's  constant.  The  usual  Weiss  molecular  field  theory  results 
from  the  assumption  that  the  local  field  has  the  Lorentz  form: 

= E + XP  , (3) 

where  E is  the  applied  field  and  X is  a constant.  The  GMFT  is  obtained1 
when  one  replaces  the  linear  Lorentz  correction  XP  by  an  internal  field 
containing  higher  powers  of  P: 

E^  = E + XP  + X'P3  + X"  P5  + ...  (4) 

1 2 

As  we  have  noted  elsewhere  ’ , this  form  for  the  local  field  is  closely 

4 

related  to  that  introduced  by  Onsager  in  treating  polar  fluids.  It  is 
readily  shown'1'  that  Eqs.  (l),  (2),  and  (4)  lead  to  a free  energy  of  the 
form 

G = l Ao(T-To)  P2  + £ Bo(T-T1)  P4  + \ CjT-T,,)  P6  + . . . , (5) 


where 
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Ao  = 3k/"/  T0  = VA0 

Bo  = 3Ao/5  (%)2  Tx  = \'/Bo  (6) 

Co  = 99  A0/175(N'J‘)4  T2  = X"/C0 

and  so  forth.  The  field  E and  susceptibility  X are  obtained  from  G by 
differentiation: 

E = 3G/SP  (7) 

and 

l/x  = dE/dP  = d2G/dP2  . (8) 

The  feature  of  the  free  energy  (5)  which  is  unusual  for  a molecular  field 

theory  is  that  the  system  may  enter  a spontaneously  polarized  state  via 

a phase  transition  which  is  either  second  order  (for  < Tq,  occurring 

at  T = T ) or  first  order  (for  T,  > T , occurring  at  T = T > T ). 
o 1 o co 

If  the  power  series  (4)  is  regarded  as  only  an  approximation  for  E^, 
it  is  reasonable  to  examine  other  functional  forms  which  yield  (4)  for 
small  values  of  P.  One  such  form,  obtained  from  the  theory  of  Pade  approxi- 
mants,  isX 

E,  = E + — — -0-  . (9) 

v l-U'/l)?1- 

Since  this  form  predicts  X'  in  good  agreement"''  with  the  value  deduced 
from  experiment  in  the  c e of  BaTiO^,  we  have  performed  the  present  calcu- 
lations using  both  forms  (4)  and  (9)  for  the  local  field. 

We  have  found  it  convenient  to  introduce  reduced  (dimensionless) 


variables 
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p = P/Np, 

e = E/N^X  (10) 

t = T/To 

(Note  that  the  unit  electric  field  Np,X  is  the  Lorentz  local  field  (3) 
at  zero  applied  field  and  saturation  polarization).  In  terns  of  the 
reduced  variables  (10),  Eqs.  (l)  and  (2)  take  the  form 


where 


p = coth  x - l/x  , 


x - 3e ./t  = 3[e+f(p)]/t 


(11) 

(12) 


The  local  field  e^  is  obtained  from  either  of  Eqs.  (4)  or  (9)  so  that 


f(p)  = 


p + ^ t.^  +.  tg  p'1  + . . . , power  series 
— P , Pade  approximant  , 

!-f  t/ 


(13) 


where  we  have  used  the  obvious  notation  t^  = T^/Tq.  If  the  free  energy 
is  normalized  a3 


8 = 


G 

3NkT 

o 


(14) 


Eq.  (5)  becomes 

g = |(t-l)  p2  + £ P4  + \ co(t-t2)  p6  + ...  , (15) 

where  bQ  = 3/5,  cq  = 99/175,  etc.  When  g is  defined  as  in  (15),  the 
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results  analogous  to  Eqs.  (7)  and  (8)  are 

e = dg/dp  (l6) 

and 

l/\X  = de/dp  = d2g/dp2  ( 17 ) 

The  pyroelectric  coefficient  9 = -(dP/dT)E  and  the  response  figure  of 
merit  for  a pyroelectric  detector  = <p/x  can  also  be  expressed  in 
terms  of  the  reduced  variables.  This  will  be  done  in  Section  IV. 
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III  STATIC  PROPERTIES 

A.  Spontaneous  Polarization, 

In  order  to  obtain  the  spontaneous  polarization  p (i.e.  the  pola- 

rization  at  zero  applied  field)  as  a function  of  temperature,  it  is 

necessary  to  solve  the  implicit  relations  (ll)-(l3)  with  e = 0.  This 

must  be  done  numerically,  and  two  simple  observations  have  proven 

most  helpful  here.  First,  for  reasons  which  will  become  clear  shortly, 

it  is  convenient  to  regard  p rather  t as  the  independent  variable^. 

Second,  Eq.  (ll)  can  then  be  solved  numerically  for  x without  regard 

to  the  form  chosen  for  e^  in  Eqs.  (12)  and  (13).  The  results  for  x as 

a function  of  p are  shown  in  Fig.  1.  In  the  numerical  calculations, 

x was  obtained  by  successive  approximation,  using  Newton's  method,  until 

Eq.  (ll)  was  satisfied  co  an  accuracy  of  one  part  in  10^.  Values  of  x, 

accurate  to  at  least  one  part  in  10  , were  obtained  for  0 < p < 1 at 

intervals  of  Ap  = 0.001.  These  results  can  be  used  to  obtain  t from 

Eqs.  (12)  and  (l3),  once  a choice  has  been  made  for  the  form  of  e,  and 

the  values  of  the  t . Figure  2 shows  p = P /Np  as  a function  of  tempe- 

rature,  using  the  power  series  form  of  e^  for  various  values  of  t^  with 

t.  = 0 for  i > 1.  As  noted  above,  when  t,  < 1 the  transition  is  second 

order,  with  pg  going  continuously  to  zero  as  t increases  toward  1.  For 

t,  > 1,  however,  there  continue  to  exist  nonzero  solutions  for  p even  for 
1 s 

t > 1.  There  are,  in  fact,  two  such  solutions,  but  it  is  easily  shown  that 

the  smaller  p , indicated  by  the  broken  curve  in  Fig.  2,  corresponds  to  a 

negative  susceptibility  (see  Section  IV)  and  so  is  not  stable.  Although 

p is  a double-valued  function  of  t,  the  inverse  relation,  t as  a function 

of  p , is  always  single-valued,  making  the  latter  problem  far  easier  to 
s 

solve  by  numerical  methods.  As  the  temperature  is  increased,  a point  is 
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reached  where  the  two  solutions  for  p coincide,  and  above  this  temperature 

S 

there  are  no  solutions.  This  is  the  transition  temperature  (Curie  tempe- 
rature) tc,  and  since  the  corresponding  spontaneous  polarization  pgc  is 

nonzero,  the  transition  produces  a discontinuous  change  in  p and  so  is 

s 

first  order. 

Essentially  similar  results  are  obtained  for  the  Pade  approximant 
form  of  e^,  as  shown  in  Fig.  3.  Since  when  p is  small  the  Pade  approximant  re- 
produces the  leading  terms  of  the  power  series  form  of  e^  with  the  same  value  of 
t , the  second-order  transitions  (t^  < l)  exhibit  essentially  identical  behavior 
in  the  two  cases  for  p <0.2.  For  larger  polarizations,  however,  the  increase 
in  p with  decreasing  temperature  for  a given  t is  always  steeper  in  the 

S X 

Pade  approximant  case  (the  only  exception  occurs  for  t^  = 0,  when  both 

forms  of  e^  are  identical  to  the  Lorentz  local  field,  and  the  conventional 

Weiss  molecular  field  results  ere  obtained).  The  two  forms  of  e^  can 

produce  dramatically  different  results  in  first-order  trc  itions  (t^  > l), 

since  only  small  increases  in  t^  beyond  1 (t^  > 1.03  in  the  Pade  case, 

t >1.08  for  the  power  series)  are  sufficient  to  make  p >0.2. 

X sc 

The  results  shown  in  Figs.  2 and  3 are  exact  in  the  sense  that  they 
were  obtained  from  numerical  so1,itions  of  the  exact  equations  (ll)-(l3). 

It  is  useful  to  compare  them  with  results  obtained  using  the  approxima- 
tions made  in  the  phenomenological  theory  of  ferroelectrxcity  introduced 

3 6 

by  Devonshire0  and  widely  used  in  fitting  experimental  data  . Devonshire 

assumed  that  the  free  energy  can  be  expanded  '3  a power  series  in  the 

7 

order  parameter  as  suggested  by  Landau  . Sir.  e the  polarization  p is  the 
order  parameter,  we  may  rite 
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X 2 x.^-L^“^_l_  /iQV 

g = 2 ap  + 4 bp  + £ cp  + . . . , (18) 

"hen  the  coefficients  a,b,c,...  are  temperature  dependent.  In  order  to 
obtain  a Curie-Weiss  law  for  the  susceptibility  (17)  above  the  Curie 
temperature,  one  must  assume  a = aQ(t-l).  Two  further  assumptions  are 
now  made  in  *'  j Devonshire  theory.  The  first  is  that  the  remaining 
coefficients  b,c,...  in  Eq.  (18)  are  temperatures  independent.  This  is 
not  the  case  in  the  GMFT,  but  including  the  appropriate  temperature 

dependence  of  these  coefficients,  as  given  by  Eq.  (15),  introduces  no 

\ 

major  complications.  The  second  assumption  (which  we  refer  to  as  the 
"Devonshire  approximation")  is  that  no  serious  error  is  made  by  neglecting 
terms  in  the  power  series  (l8)  of  order  higher  than  p^.  It  is  the  conse- 
quences of  this  second  assumption  which  we  wish  to  examine.  When  the 
series  is  truncated,  Eq.  (16)  is  used  to  obtain  e,  and  e is  then  set 
equal  to  zero,  we  obtain  for  the  spontaneous  polarization 


- - h (1  * i/^) 


(19) 


For  thu  GMPT,  Eq.  (19)  becomes 


2 356i  r 

ps  - - ssr  a ±i1- 


44  6o62, 

7 ,2*  ’ 


(20) 


where  6^  = t-t^  and  t = 1.  For  t^  < 1 only  one  choice  of  sign  in  (20)  gives 

2 

a real,  positive  value  of  p . For  t > 1 (first-order  transitions)  one 

s X 

p 

obtains,  as  in  the  exact  solution,  only  one  real,  positive  value  of  pg  for 

p 

0 < t < 1,  but  two  such  values  for  1 < t < tQ.  Again,  the  smaller  pg  leads  . 

to  a negative  susceptibility  and  is  unphysical.  Results  for  p as  a function 

s 
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of  t with  tg  = 0 for  several  values  of  t^  are  shown  in  Fig.  4,  where  the 

solid  curves  were  calculated  using  the  Devonshire  approximation  and  the 

exact  results  (broken  curves)  are  given  for  comparison.  One  can  readily 

see  that,  although  the  results  of  the  Devonshire  approximation  approach 

the  exact  solution  for  small  p , significant  differences  can  appear  for 

s 

p >0.4,  and  the  discrepancy  at  a given  value  of  t increases  as  t 
S 

increases.  When  t^  is  greater  than  one,  so  that  the  phase  transition  is 
first  order,  the  polarization  near  the  critical  temperature  quickly 
becomes  so  large  that  the  results  of  the  Devonshire  approximation  are  no 
longer  valid  even  close  to  the  transition  temperature  (note  the  curves 
for  t^=2  in  Fig.  4).  Although  the  present  results  apply  only  to  the  GMFT, 
they  strongly  suggest  that  the  Devonshire  approximation  should  be  used 
with  caution  when  dealing  with  first-order  phase  transitions  and  that 
significant  errors  may  occur  in  quantities  calculated  from  an  expansion 
of  the  free  energy  which  has  been  truncated  in  the  usual  way. 

B.  Critical  Temperature  and  Polarization. 

Both  t = T /T  and  p = P /Nu,  can  be  obtained  by  locating  the  maximum 
c c'  o sc  sc' 

in  t as  a function  of  p„,  and  both  increase  with  increasing  t as  shown  in 

£ -L 

Figs.  5 and  6.  A more  accurate  method  of  determining  t and  p has  been 

c sc 

found,  however,  and  was  used  to  calculate  the  values  shown  in  the  figures. 
Tills  method  is  based  on  results  for  the  pyroelectric  coefficient  and  so  will 
be  discussed  in  Section  IV. 

An  important  point  to  be  noted  here  is  that  the  singularity  in  the 
Pade  approximant  form  of  e^  for  t^  > 5/3  causes  tc  to  diverge  and  pgc  to 
approach  1 as  t^  increases  toward  5/3  (indicated  by  the  broken  vertical 
line  in  Fig.  6).  Hence,  in  this  case,  only  values  of  t < 5/3  can  be 
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regarded  as  physically  significant.  It  is  interesting  to  note  that  when 


the 


power  series  form  of  e^  is  used  p^^  never  reaches  1,  as  it  does  for 
the  Pad^  approximant  form,  but  rather  approaches  a limit  of  approximately 
0. 746  as  t -*  00 , 
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IV.  DYNAMIC  jflOFERTIES 
A.  Susceptibility. 

The  reduced  inverse  susceptibility  l/\X  defined  in  Eq.  (17)  can  be 
obtained  by  differentiating  Eq.  (ll)'  with  respect  to  e and  solving  for 
de/dp.  The  result  is 


Here 


l/Xx  = 3P"  “ f'(p)  • 


p'  = dp/dx  = 1 - coth2x  + 1/x2  , 


(21) 


(22) 


where  x is  defined  in  Eq.  (12),  and 


f'(p)  = df/dp  = 


r q 2 QQ  4 

1 + | tlP  + H tgP  + . . . , power  series 
(l  + ^ t1p2)/(l  - ^ t^p2)2,  Pade  approximant , 


(23) 


where  f(p)  is  defined  in  Eq.  (13).  Equations  (2l)-(23)  have  been  used  to 
calculate  l/Xx  with  e = 0 for  various  values  of  t^,  and  the  results  are 
shown  in  Figs.  7 and  8,  for  the  power  series  and  Pade  approximant  forms  of 
e^,  respectively.  (Note  the  logarithmic  scale  for  l/Xx).  In  the  Devonshire 
approximation  we  have 


1/XX  = a + 3bp2  + 5cp4  , 


(24) 


and  for  the  GMFT  this  becomes 


(25) 


VXX  = 60  + | 6iP2  + |f  62P4  » 

where  6 = t-ti  as  in  Eq.  (20).  Using  the  Devonshire  approximation  (25), 


we  obtain  the  results  shown  in  Fig.  9.  As  one  might  expect,  all  three 
calculations  give  similar  results  when  the  polarization  is  small,  i.e.  near 
the  transition  temperature  for  second-order  transitions.  As  in  the  case  of 
the  spontaneous  polarization,  the  differences  increase  with  increasing 
polarization  and  are  larger  the  greater  the  absolute  value  of  t . It  is 
interesting  to  note  that  the  results  for  t = -1  and  t = 0 clearly  approach 
one  another  near  the  transition  temperature  in  all  three  calculations,  while 
the  results  for  t^  = 1 are  distinctly  different  in  the  temperature  range 
shown  in  Figs.  7-9-  The  reason  for  this  will  become  clear  as  we  examine 
in  detail  the  behavior  of  l/Xx  in  the  vicinity  of  the  transition  temperature. 

It  is  readily  shown  that  the  exact  result,  Eq.  (2l),  leads  to  a Curie- 
Weiss  law  in  the  paraelectric  phase.  If  all  the  quantities  in  (2l)  are 
expanded  in  power  series  in  the  polarization,  we  obtain 

1/AX  = t - 1 + 2 (t-tl)  p + ...  . (26) 

In  the  paraelectric  phase  the  polarization  p vanishes,  and  we  obtain  the 
Curie-Weiss  law, 

1/XX  = t - 1 . (27) 

If  an  electric  field  is  applied  to  the  paraelectric  phase  in  order  to  induce 
a nonzero  polarization,  Eq.  (24)  shows  that  the  change  in  l/\X  will  be  pro- 
portional to  p“,  in  agreement  with  the  Devonshire  theory''.  In  addition,  the 
GMFT  predicts’'"  that  the  coefficient  of  p^  will  be  linear  in  t and  vanish 
at  the  characteristic  temperature  t . This  behavior  has  been  observed  ex- 

g 

perimen tally  by  Drougard  el  al  and  provides  a means  for  determining  t . 
Measurements  are  now  in  progress  to  confirm  the  results  of  Ref.  8 for 
BaTiO^  (t^  2 1-17)  and  to  obtain  values  of  t^  for  other  materials  of  interest. 
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In  materials  which  undergo  second-order  transitions,  p goes  conti- 
nuously to  zero  as  t -*  1 from  below,  so  that  close  to  the  transition 
temperature  p can  be  arbitrarily  small.  Hence  there  exists  a range  of 
temperature  in  which  it  is  possible  to  assume  a more  drastic  form  of  the 

Devonshire  approximation,  namely  that  the  term  in  the  free  energy  (l8) 

6 3 

proportional  to  p can  be  neglected  as  well  as  all  the  higher-order  terms  . 

The  sDontaneous  polarization  is  then  given  by  p = - a/b,  which  for  the 
- - s 

GMFT  becomes 


2 

P„  = " 


(28) 


For  e = o we  can  substitute  (?8)  in  (26)  and  obtain 


n / *\  * . _ nf  4-  n N 

**“•/  * */v  ~ — \ — / f 


(0Q\ 


■3 

which  is  the  well-known  result^  that  in  second-order  transitions  the 
susceptibility  has  a Curie-Weiss  behavior  below  as  well  as  above  the  tran- 
sition temperature,  with  the  ratio  of  the  slopes  equal  to  -2,  Figure  10 
snows  the  results  of  the  exact  caxmuafion  using  Eqs.  (2l)-(23)  for  se- 
veral second-order  transitions  (^  < !)•  0nJ-y  the  region  near  the  tran- 
sition temperature  is  shown,  and  the  scale  for  l/Xx  is  now  linear.  In 
this  temperature  range,  the  polarization  is  small  enough  that  the  power 
Series  anu  Fade  approximant  forms  of  e^,  as  well  0.G  ulic  i/cVGiibii ii'c  appro ^ 
ximaticn,  yield  essentially  identical  results  except  for  the  case  of 
t^  = 1 which  will  be  discussed  below.  Since  the  result  (27)  is  exact  for 
the  paraelectric  phase,  the  behavior  of  l/Xx  is  the  same  for  both  forms 
of  p.  and  an  values  of  t. . In  the  ferroelectric  phase,  Eq,  (29)  is 

i - ■ 

2 

correct  through  terms  of  order  p , so  that  its  validity  can  extend  over 
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a substantial  temperature  range,  especially  for  small  values  of  tx  which, 
as  noted  above,  lead  to  a relatively  gradual  Increase  in  polarization  as 

the  temperature  is  lowered. 

The  case  ^ = 1 must  be  treated  separately  because  Its  critical  behavior 
is  different  from  that  of  other  second-order  transitions.  Equation  (26),  which 
predicts  that  the  spontaneous  polarization  has  the  temperature  dependence 
D ^ « i_t,  fails  in  this  case.  Using  Eq.  (20)  with  ~ \ _ 1 one 

* g ' | 

can  show  that  the  critical  behavior  is  actually  Ps  * 1-t.  One  effect  of 
this  change  in  critical  exponent  is  readily  seen  in  Figs.  2-4,  where  the 
curvature  of  t as  a function  of  Ps  vanishes  at  ^ = 1.  The  unique  critical 
behavior  of  l/XX  for  t^l  is  clearly  visible  in  Fig.  10,  where  the  Fade' 
approximant  form  of  e^  (solid  curve),  power  series  form  of  e^  (broken  curve), 
and  Devonshire  approximation  (dash-dot  curve),  yield  significantly  differen 
results.  Unusual  critical  behavior  for  the  case  tfl  should  not  be  sur- 
prising, since  the  point  (p.t,^)  = (0,1, l)  is  a tricritical  point  in  the 
sense  of  Griffiths9,  and  the  critical  exponents  at  such  a point  may  differ 
from  those  at  neighboring  points.  A detailed  study  of  the  critical  beh 

for  t =1  is  now  in  progress. 

A final  word  should  be  said  about  the  low-temperature  behavior  of 
l/XX.  In  the  exact  results  (Fig.  7 and  8),  l/XX  diverges  as  t - 0,  re- 
flecting the  fact  that  X - 0 as  the  polarization  saturates.  For  the 
Devonshire  approximation,  the  results  are  rather  peculiar  and  unphysical. 

It  can  be  shown  with  the  aid  of  Eq.  (20)  that  for  < 0 and  tg  = 0,  we 

have  p 2 - -5/3t  as  t - 0.  Thus  P0  saturates,  but  at  a value  which  depends 
on  tx  rather  than  at  ?s  = 1.  If  this  finite  x^it  for  ?s  is  substituted 
in  Eq.  (24),  one  obtains  a finite  limit  for  l/XX  rather  than  a divergence 
(see  the  curve  for  t±  = -1  in  Fig.  9).  On  the  other  hand,  if  ^ > 0, 
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both  the  polarization  (20)  and  the  inverse  susceptibility  (24)  diverge 
as  t -»  0.  We  conclude,  therefore,  that  the  Devonshire  approximation  must 
always  fail  at  low  temperatures:  either  the  polarization  or  the  suscepti- 
bility will  behave  in  an  unphysical  manner. 

B.  iyroelectric  Coefficient. 

In  terms  of  the  reduced  variables  (10),  one  can  define  a reduced 
pyroelectric  coefficient 

cp*  = - (9p/3t)e  , (30) 

which  is  related  to  the  conventional  pyroelectric  coefficient 
cp  = - (3P/3T)E  by 

(3D 

* 

In  the  following  we  discuss  cp  only  for  e = 0,  which  corresponds  to  the 

usual  experimental  situation.  It  should  be  noted,  however,  that  an  applied 

* 

electric  field  produces  significant  changes  in  cp  and  can  induce  a pyro- 
electric effect  above  the  Curie  temperature^0.  These  effects  will  be 
discussed  in  a future  report. 

. * . . 

An  exact  expression  for  cp  can  be  obtained  by  differentiating  Eq.(ll) 

wJth  respect  to  t and  solving  for  dp/dt.  The  result  is 


* 

cp 


xp  ’ 

t-3p'f”(p) 


(32) 


where  p'  and  f'  are  defined  in  Eqs.  (22)  and  (23).  The  pyroelectric 

coefficient  calculated  from  (32)  with  e = 0 (p=p  ) is  shown  as  a function  of 

s 
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reduced  temperature  in  Figs.  11  and  12  for  various  values  of  the  charac- 
teristic temperature  t^,  using  the  power  series  and  Pade/  approximant  forms 
of  the  local  field,  respectively. 

In  the  Devonshire  approximation,  the  reduced  pyroelectric  coefficient 
for  e = 0 is1 


2 i* 

„ a'+b'p  +c'p 
cp  - § JL 

2p  (b+2cp 

P o 


(33) 


where  the  prime  indicates  differentiation  with  respect  to  t.  For  the  GMFT 
this  becomes 


* -K 

9 = 


m 


2h  (~  6 
* s ' b 1 


3 A + M 6 p 2) 
17b  2 ps  ' 


(34) 


where  6^^  = t-t  as  in  Eq.  (20).  The  results  of  Eq.  (34)  are  shown  in 
Fig.  13  for  the  same  values  of  t used  in  the  calculations  of  Fig.  11  for 
the  power  series  form  of  e^. 

The  exact  results  shown  in  Figs.  11  and  12  are  fully  consistent 

with  the  polarization  curves  in  Figs.  2 and  3-  Large  valuer  uf  t produce 

* . 

a flat  polarization  curve  (and  hence  small  cp  ) at  low  temperatures,  followed 

* 

by  a sharp  drop  in  pg  giving  a large  cp  near  the  transition  temperature. 

For  small  t , cp  is  relatively  larger  at  low  temperature  and  smaller  near 

-x 

the  transition,  giving  rise  to  the  crossings  of  the  cp  curves  for  t < 1 

which  characterize  the  results  in  Figs.  11  and  12.  Note  also  that  the 

steeper  polarization  curves  produced  by  the  Fade  approximant  form  of  e^ 

* 

lead  to  smaller  values  of  cp  at  low  temperatures  and  larger  values  near 


the  transition  than  are  obtained  with  the  power  series  form. 
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Here  again  the  Devonshire  approximation  fails  to  give  results  con- 
sistent with  the  exact  solutions  (see  Fig.  13).  For  t^  > 0,  the  pyro- 
electric coefficient  diverges  as  t -*  0.  This  is  due  to  the  unphysical 

divergence  of  p noted  above  in  connection  with  the  susceptibility.  Un- 
s 

fortunately,  this  low- temperature  divergence  is  stronger  than  the  one  at 

* 

the  transition  temperature  and  dominates  the  behavior  of  cp  throughout 
most  of  the  temperature  range.  The  crossings  obtained  in  the  exact  results 
are  eliminated,  and  the  range  in  which  the  approximation  yields  satisfac- 
tory results  is  even  smaller  than  in  the  case  of  the  polarization  or 
susceptibility. 

Equation  (32)  for  the  pyroelectric  coefficient  provides  the  basis 

for  an  accurate  calculation  of  the  temperature  t at  which  a first-order 

transition  occurs  when  t > 1.  As  can  be  seen  in  Figs.  2-4,  tR  corresponds 

to  the  ''nose"  of  the  polarization  curve,  beyond  which  no  nonzero  solutions 

of  Eqs.  (ll)-(l3)  exist.  At  this  point  dp  /dt  diverges.  This  cannot  happen 

s 

for  p <1  unless  the  denominator  of  Eq.  (32)  vanishes.  Thus  the  transition 
s 

temperature  tQ  and  the  spontaneous  polarization  at  that  temperature  psc  are 
related  by 


t 

c 


3 p'  f'(p  ) , 

■^sc  Visc  ’ 


(35) 


where  p'  is  simply  p'  as  defined  in  Eq.  (22)  evaluated  at  e = 0 and 
s c 

t = tc,  and  f'  is  defined  in  Eq.  (23).  Ifow  Eq,(l2)  gives  a second  relation 

between  t and  p 
c sc 


t 

c 


3 f(p  )/x  , 

v±sc//  sc  * 


(36) 


where,  as  usual,  the  subscripts  s and  c indicate  e = 0 and  t = tc,  respec- 
tively. We  have  found  it  most  convenient  to  eliminate  t from  Eqs.  (35)  and 


T-*  1 m >>■■  11 
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(36)  and  to  solve  for  as  a function  of  p (assuming,  in  the  power  series 

X sc 

form  of  e,,  that  t,  = 0 for  i > l).  This  gives 

■C  a. 


where 


■ sc 


- X p' 

sc  sc 


r p 


sc 


(37) 


3 x p'  - p , power  series 
sc  sc  sc 


D = 


x p'  + p , Pade  approximant  . 
SC  sc  sc 


(38) 


By  regarding  p as  the  independent  variable,  solving  (37)  and  (38)  for 
s c 

t and  then  evaluating  tc  from  (36),  one  obtains  the  results  shown  in 
Figs.  5 and  6. 


C.  Response  Figure  of  Merit 

It  is  well-known  that  the  sensitivity  of  a pyroelectric  detector  is 
proportional  to  the  ratio  M1  = <p/x,  often  referred  to  as  the  response  figure 
of  merit.  In  terms  of  reduced  variables,  the  corresponding  ratio  is 


which  is  related  to  the  conventional  figure  of  merit  by 


m 


1 


Combining  (2l)  and  (32)  gives  the  simple  result, 


(39) 


(4o) 


n»1  = x/3.  (4l) 

In  tne  following  we  will  discuss  the  figure  of  merit  at  zero  applied  electric 
field.  This  is  the  simplest  case  and  corresponds  to  the  situation  encoun- 
tered in  most  detector  applications,  but  it  should  be  noted  that  pyroelectric 
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detectors  have  been  constructed  which  require  a d.c.  bias  field,  and  the 

characteristics  of  this  mode  r*  operation  have  been  analyzed . ^ 

Using  Eq.  (4l),  we  have  calculated  m^  as  a function  of  temperature 

for  various  values  of  t^,  using  the  power  series  form  of  e^  (with 

t^  = 0 for  i > l)  and  using  the  Pad!  approximant  form.  The  results,  shown 

in  Figs.  14  and  15,  are  quite  similar  for  the  two  forms  of  the  local  field. 

The  divergence  of  as  t H 0 is  due  to  the  fact  that  cp*  remains  finite 

while  x goes  to  zero.  This  results  from  the  assumption,  implicit  in 

Eq.  (l),  that  the  elementary  dipoles  can  assume  any  orientation  in  space, 

that  is  they  behave  "classically"  as  in  the  "infinite-spin"  limit  of  the 

Weiss  theory  of  ferromagnetism.  Other,  equally  plausible,  assumptions  lead 

to  the  result  that  <p*  -*  0 as  t -*  0.  We  will  return  to  this  point  shortly. 

From  a practical  point  of  view,  this  enhancement  of  may  not  be  useful, 

even  if  correct,  simply  because  the  total  susceptibility  can  never  vanish, 

but  will  instead  be  dominated,  as  p saturates,  by  the  electronic  and  non- 

s 

ferroelectric  lattice  contributions  which  are  neglected  in  the  present  model. 
In  the  Devonshire  approximation,  Eq.  (39)  becomes'*' 


m±  = a'ps  + b'pg3  + c'ps5  , (42) 

where  the  prime  again  indicates  differentiation  with  respect  to  t.  For 
the  GMFT,  Eq.  (42)  yields 


m. 


<1+ 1 


+ 


22  P h) 

175  s ' 


9 


(43) 


from  wi-ich  we  have  calculated  the  results  shown  in  Fig.  16,  For  small 
values  of  t^  the  Devonshire  approximation  is  in  reasonable  agreement  with 
the  exact  results  at  sufficiently  high  temperatures,  indicating  considerable 
cancellation  of  the  errors  occurring  in  cp*  and  l/Xx  separately.  This  might 


-21 


be  expected  from  the  fact  that  Eq.  (43)  is  correct  up  to  terms  of  order  p . 

s 

There  remains  one  extremely  important  point  to  be  made  concerning 

Eq.  (4l) . Since  depends  only  on  x,  which  in  turn  depends  only  on  pg 

via  Eq.  (ll),  we  conclude  that  m depends  only  on  p . In  this  case,  re- 

garding  the  polarization  as  the  independent  variable  is  more  than  the 

computational  convenience  it  has  been  up  to  this  point.  Regarding  m^  as 

a function  of  p predicts  a universal  behavior  of  the  figure  of  merit  as 
s 

a function  of  the  spontaneous  polarization,  independent  of  the  choice  of 

a form  for  the  local  field  (and  so  necessarily  independent  of  the  values 

of  the  t^).  This  is  shown  in  Fig.  17,  where  m^  is  given  as  a function  of 

1-p  for  the  "classical"  case  described  by  p = coth  x - l/x  (see  curve 
s 

labelled  J = “>).  The  variable  1-p  , which  increases  with  temperature,  is 

s 

used  rather  than  p so  that  comparison  can  be  made  with  Figs.  l4-l6.  In 
s 

-udition,  we  have  calculated  m^  for  the  "spin-^"  case  in  which  the  pola- 
rization is  given  by  p = tanh  x,  and  for  the  Devonshire  approximation  to 
the  "classical"  case.  These  are  indicated  in  Fig.  17  by  the  curve  labelled 
J = \ and  the  broken  curve,  respectively. (Note  that  for  J = ^ we  find  that 
m^  remains  finite  as  t -*  0.  This  is  because  cp*  -*  0 as  t -*  0 in  this  case.) 
Perhaps  surprisingly,  the  three  results  are  indistinguishable  for  p <0.5. 

S *■“ 

So  while  Eq.  (4l)  provides  a direct  test  of  whether  the  dipoles  behave 
classically  or  have  only  certain  allowed  orientations,  this  information 
can  be  obtained  only  from  measurements  made  near  saturation  polarization, 
that  is  in  the  low-temperature  region  which  has  not  previously  been  of 
sufficient  interest  to  warrant  detailed  experimental  study.  Low-temperature 
measurements  of  the  pyroelectric  coefficient  and  susceptibility  (corrected 
for  electronic  and  non-ferroelectric  lattice  contributions)  are  planned  for 


the  near  future. 
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Equations  ( 4o)  and  (4l)  predict  that  for  a given  polarization  the 
experimental  response  figure  of  merit  should  scale  from  one  material 
to  another  as  Np\/To.  Data  is  now  being  collected  to  check  this  prediction 
which,  if  verified,  could  serve  as  a useful  guide  in  obtaining  larger  values 
of  M . 

Finally,  it  is  interesting  to  note  that  the  scale  factor  in  Eq.  (4o) 
can  be  written,  with  the  aid  of  Eq.  (6),  as 


T 

0 _ JJ;_ 

NpX  = 3k  ' 


(44) 


Thus  the  experimental  figure  of  merit  is  inversely  proportional  to  the 
dipole  moment  p.  Since  the  left-hand  side  of  Eq.  (44)  can  be  calculated 
from  experimental  quantities  (saturation  polarization,  Curie  coefficient, 
and  Curie-Weiss  temperature),  we  can  now  determine  p indirectly.  The  value 
of  Eq.  (44)  may  be  greatly  increased  if  it  proves  possible  to  calculate  p 
directly  from  lattice  dynamical  models. 
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V.  COMPARISON  WITH  EXPERIMENT 

A detailed  comparison  of  the  predictions  of  the  GMFT  presented  above 
with  the  experimental  data  available  in  the  literature  is  now  in  progress 
and  will  be  the  subject  of  a future  report.  The  preliminary  phase  of  this 
comparison  has  already  yielded  some  important  information  and  conclusions 
which  will  be  described  here. 

As  noted  in  Section  IVA,  the  characteristic  temperature  t^  can  be 

determined  from  the  electric  field  dependence  of  l/X-X.  This  has  been  done 

to  date  for  only  one  material,  however,  so  that  it  would  be  most  valuable 

to  have  another  method  of  finding  t^  from  data  already  available  in  the 

literature.  In  the  case  of  a first-order  phase  transition,  the  transition 

temperature  tc  depends  on  t^  as  shown  in  Figs.  5 and  j for  the  power  series 

and  Pade  approximant  forms  of  the  local  field,  respectively.  Hence,  in  this 

case  t can  be  determined  from  the  observed  value  of  tc>  Unfortunately, 

many  materials  of  interest,  such  as  triglycine  sulfate  (TGS)  and  its  iso- 

morphs,  undergo  second-order  transitions  which  occur  at  t = 1 independent 

of  the  value  of  t^.  So  yet  another  method  of  finding  t^  is  needed,  and  our 

preliminary  results  indicate  that  fitting  the  temperature  dependence  of  the 

spontaneous  polarization  p„  can  be  used  for  this  purpose. 

s 

12 

Figure  It)  shows  the  experimental  results  of  Triebwasser  for  the 
spontaneous  polarization  of  TGS  as  a function  of  temperature  (both  in 
reduced  units).  'The  solid  curve  is  the  prediction  of  the  GMFT,  using  the 
Pade  pproxiinant  form  of  the  local  field  with  t^  = O.85.  The  agreement 
between  theory  and  experiment  is  well  within  the  experimental  error.  The 
value  t = 0 85  was  determined  by  trial  and  error,  using  variations  of 
At^  = 0.05.  The  fit  for  t^  = O.85  is  significant^  oetter  than  for 
tj  = O.JC  or  t,  = 0.90,  but  the  amount  of  computation  required  for  further 
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refinement  is  too  great  to  be  performed  by  hand.  A computer  program  is  now 
being  written  to  carry  out  a least-squares  fit  and  determine  t to  the 
greatest  precision  compatible  with  the  accuracy  of  the  experimental  data. 

Our  preliminary  result,  t = 0.85  ± 0.03  for  TGS,  can  also  be  obtained 
from  considerations  based  on  the  Devonshire  approximation.  The  more  drastic 

form  of  this  approximation,  neglecting  terms  in  the  free  energy  (l8)  of  order 

c 

p and  higher,  leads  (for  the  GMFT)  to  Eq.  (28)  which  can  be  solved  for  t 
a function  of  p and  t: 

t = t + -2-g  (t-1)  (45) 

3ps 

Averaging  the  three  experimental  points  in  Fig.  18  having  the  smallest 
values  of  p yields  t =0.84,  in  good  agreement  with  the  trial-and-error 
fit  shown  in  Fig.  18. 

Even  though  the  computer  program  required  to  obtain  optimum  values 

of  t from  the  observed  temperature  dependence  of  p is  not  yet  completed, 

X s 

it  is  already  possible  to  draw  certain  conclusions  from  the  preliminary 

fit  for  TGS.  First,  the  Pade  approximant  form  of  the  local  field  e^  gives 

a significantly  better  fit  than  does  the  power  series  form.  The  value 

t ^0-8^  obtained  for  small  p from  the  Devonshire  approximation  is  inde- 

pendent  of  the  form  chosen  for  e^,  but  when  this  value  is  used  to  calculate 

p for  small  t (i.e.  large  p ),  the  Pade  approximant  form  of  e.  gives 
s s 

results  in  good  agreement  with  experiment  (see  Fig.  18)  while  the  power  series 

form  with  t^  =0  for  i > 1 results  in  large  deviations  from  experiment  when 

p is  greater  than  approximately  0.5. 
s 

Second,  the  relatively  large  value  of  t^  obtained  for  TGS  confirms 
the  superiority  of  the  GMFT  to  the  conventional  molecular  field  theory 
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(which  is  identical  to  the  GMFT  with  t^  = 0 for  all  i).  In  the  context  of 
the  GMFT,  the  phase  transition  in  TGS  is  "nearly  first  order"  since  t^  > 1 
results  in  first-order  transitions.  The  onset  of  the  spontaneous  polariza- 
tion for  t < 1 is  fhr  steeper  than  predicted  by  the  conventional  molecular 
field  theoiy  and  cannot  be  accounted  for  without  introducing  the  gene- 
ralized molecular  field  described  by  Eq.  (4'. 

Finally,  the  success  of  the  GMFT  in  fitting  the  spontaneous  polari- 
zation of  TGS  encourages  us  to  believe  that  the  susceptibility,  pyroelectric 
coefficient,  and  figure  of  merit  of  this  material,  as  well  as  the  properties 
of  other  ferroelectrics,  can  be  predicted  from  the  GMFT.  An  extensive 
program,  mentioned  above,  is  now  in  progress  to  compare  the  predictions 
of  the  GMFT  with  published  expei imental  results  and  with  those  obtained 


in  our  laboratory. 
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VI.  SUMMABY 

The  investigations  described  in  this  report  were  undertaken  with  the 

goal  of  elucidating  the  physical  mechanism  responsible  for  the  pyroelectric 

effect  in  ferroelectric  crystals,  in  the  hope  that  an  understanding  of 

this  mechanism  would  lead  to  practical  guidelines  for  the  selection  of 

infrared  detector  materials  and  to  greater  insight  into  the  factors  which 

limit  the  performance  of  pyroelectric  detectors.  We  have  reported  pre- 
1 2 

viously  * on  an  important  generalization  of  the  well-known  molecular 
field  theory  of  ferroelectricity.  (These  results  are  summarized  briefly 
in  Section  II.)  With  this  generalization  it  is  possible  to  overcome  the 
major  objections  which  have  been  raised  to  the  use  of  a molecular  field 
theory  in  this  context^.  Much  of  the  work  reported  here  is  of  an  essen- 
tially technical  nature:  the  generalized  molecvlar  field  theory  (GMFT) 
is  used  to  obtain  exact  expressions  for  the  ferroelectric  properties  of 
crystals  (spontaneous  polarization,  dielectric  susceptibility,  critical 
temperature  and  polarization)  and  for  the  resulting  pyroelectric  properties 
(pyroelectric  coefficient  and  response  figure  of  merit  for  radiation 
detection).  These  results  are  presented  in  Sections  III  and  IV  for  various 
values  of  the  parameters  which  appear  in  the  GMFT.  An  extensive  comparison 
of  these  predictions  with  the  available  experimental  data  is  now  under  way, 
and  some  highly  encouraging  results  obtained  in  the  preliminary  stages  of 
this  work  are  reported  in  Section  V. 

Two  major  qualitative  results  have  emerged,  however,  from  these  quan- 
titative considerations.  In  Section  II  we  introduce  a set  of  reduced 
(dimensionless)  variables  which  greatly  simplify  the  calculations,  and  in 
Section  III  we  find  that  regarding  the  polarization  rather  than  the  tempe- 
rature as  the  independent  variable  results  in  a large  saving  of  computational 
effort.  When  we  turn  our  attention  to  the  response  figure  of  merit  for 
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pyroelectric  detectors  (Section  IV) , we  find  that  these  two  devices  are 
not  merely  convenient  computational  aids.  Rather,  they  lead  directly  to 
a striking  universal  behavior  for  the  figure  of  merit.  We  find  that  the 
figure  of  merit,  in  reduced  units,  is  always  given  by  the  same  universal 
function  of  the  polarization,  independent  of  the  many  details  of  the  GMFT 
which  can  vary  from  one  material  to  another  and  affect  the  behavior  of 
other  properties  such  as  the  susceptibility  or  pyroelectric  coefficient. 
Work  is  now  in  progress  to  verify  experimentally  this  universal  behavior 
and  to  determine  whether  it  can  be  exploited  as  a tool  for  optimizing 
the  response  figure  of  merit. 

Also,  we  have  compared  the  exact  results  of  the  GMFT  with  the 
Devonshire  formalism  obtained  from  the  GMFT  by  truncating  the  power-series 
expansion  for  the  free  energy.  Our  results  show  that  in  general  the 
Devonshire  formalism  is  a poor  approximation  to  the  GMFT,  especially  for 
first-order  transitions.  Only  in  the  case  of  the  response  figure  of  merit 
are  the  results  of  this  approximation  accurate  over  a wide  enough  range 
of  polarization  to  be  considered  useful. 
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Figure  Captions 

Fig. 1 - Solutions  of  the  polarization  equation  p = coth  x - l/x. 

Fig. 2 - The  spontaneous  polarization  p = P /Np,  as  a function  of  the  reduced 

s s 

temperature  t = T/Tq  for  the  power  series  form  of  the  local  field 
e^  and  various  values  of  the  characteristic  temperature  t^  = T^/Tq 
(with  t = 0 for  i > l). 

Fig. 3 - The  spontaneous  polarization  p = P /Nja  as  a function  of  the  reduced 

temperature  t = T/Tq  for  the  Pade  approximant  form  of  the  local  field 

e^  and  various  values  of  the  characteristic  temperature  t^  = T^/Tq. 

Fig. 4 - The  spontaneous  polarization  p = P /Njj,  as  a function  of  the  reduced 

s s 

temperature  t = T/Tq  calculated  in  the  Devonshire  approximation 
(solid  curves)  is  compared  with  the  exact  results  using  the  power 
series  form  of  the  local  field  e^.  The  same  values  of  the  characte- 
ristic temperature  t.^  = T^/T  are  used  as  in  Fig.  2,  with  ^ = 0 
for  i > 1. 

Fig. 5 - The  Curie  temperature  tc  = Tc/Tq  and  the  spontaneous  polarization  at 

the  Curie  temperature  p = P /N(jl  as  a function  of  the  characteristic 

SC  sc 

temperature  t^  = T^/T^  (with  t^  = 0 for  i > l)  using  the  power  series 
form  of  the  local  field  e^. 

Fig. 6 - The  Curie  temperature  t^  = Tc/Tq  and  the  spontaneous  polarization  at 

the  Curie  temperature  p = P /N^i  as  a function  of  the  characteristic 

sc  sc 

temperature  t^  = T^/T  using  the  Pade  approximant  form  of  the  local 
field  e^. 

Fig. 7 - The  reduced  inverse  susceptibility  l/\x  = de/dp  at  e = 0 (note 

logarithmic  scale)  as  a function  of  the  reduced  temperature  t = T/T 

' o 

for  the  power  series  form  of  the  local  field  e^  and  various  values  of 
the  characteristic  temperature  t^  = T^/Tq  (with  t^  = 0 for  i > l). 
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Fig. 8 - The  reduced  inverse  susceptibility  l/\X  = de/dp  at  e = 0 (note 

logarithmic  scale)  as  a function  of  the  reduced  temperature  t = T/T 

o 

for  the  Pade/  approximant  form  of  the  local  field  e^  and  various 

values  of  the  characteristic  temperature  t = T^/T  . 

Fig. 9 - The  reduced  inverse  susceptibility  l/\\  = de/dp  at  e = 0 (note 

logarithmic  scale)  as  a function  of  the  reduced  temperature  t = T/T  , 

calculated  in  the  Devonshire  approximation.  The  seme  values  of  the 

characteristic  temperature  t^  = T^/Tq  are  used  as  in  Fig.  7,  with 

t.  = 0 for  i > 1. 
l 

Fig. 10-  The  reduced  inverse  susceptibility  l/\x  = de/dp  at  e = 0 (note 

linear  scale)  as  a function  of  the  reduced  temperature  t = T/Tq 

in  the  region  near  the  transition  temperature,  for  selected  second- 

order  transitions  (t  = T^/Tq  < l).  For  t > 1 and  for  t < 1 with 

t1  < 1,  the  power  series  and  Pade  approximant  forms  of  the  local 

field  e^  give  essentially  identical  results.  For  t^  = 1 the  results 

of  the  Pad^  approximant,  power  series,  and  Devonshire  approximation 

are  indicated  by  the  solid,  broken,  and  dash-dot  curves,  respectively. 

Fig. 11-  The  reduced  pyroelectric  coefficient  eg*  = -(dp  /dt)  as  a function  of 

s 

the  reduced  temperature  t = T/Tq  for  the  power  series  form  of  the 

local  field  e^  and  various  values  of  the  characteristic  temperature 

tn  = T../T  (with  t.  = 0 for  i > l). 

1 1'  o x 

Fig. 12-  The  reduced  pyroelectric  coefficient  eg*  = -(dp  /dt)  as  a function 

s 

of  the  reduced  temperature  t = T/Tq  for  the  Pad^  approximant  form 
of  the  local  field  e^  and  various  values  of  the  character:} st-i" 
temperature  t^  = T^/To. 

Fig. 13-  The  reduced  pyroelectric  coefficient  eg*  = -(dp  /dt)  as  a function 

s 

of  the  reduced  temperature  t = T/Tq,  calculated  on  the  Devonshire 
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approximation.  The  same  valves  of  the  characteristic  temperature 

t = T /T  are  used  as  in  Fig.  11,  with  t.  = 0 for  i > 1. 

1 1'  o 

Fig. i4  - The  reduced  response  figure  of  merit  m^  = cp^'/lX  at  zero  applied 
field  as  a function  of  the  reduced  temperature  t = T/Tq  for  the 
power  series  form  of  the  local  field  e,  and  various  values  of  the 
characteristic  temperature  t^  ~ T^/T^  (with  t^  = 0 for  i > l) . 

Fig. 15  - The  reduced  response  figure  of  merit  = cp*/lX  at  zero  aPPlied 
field  as  a function  of  the  reduced  temperature  t = t/Tq  for  the 
Pade  approximant  form  of  the  local  field  e^  and  various  values  of 
the  characteristic  temperature  ^ = T^/T^. 

Fig. 16  - The  reduced  response  figure  of  merit  = cp*/lx  at  zero  applied 

field  as  a function  of  the  reduced  temperature  t = T/T  , calcula- 
ted in  the  Devonshire  approximation.  The  same  values  of  the 
characteristic  temperature  = T^T^  are  used  as  in  Fig.  14, 
with  t^  = 0 for  i > 1. 

Fig.  17  - The  reduced  response  figure  of  merit  m^  = Vx'/lX  at  zero  applied 

field  as  a function  of  1 - po , showing  its  universal  behavior.  The 
curves  J = | and  J = 00  correspond  to  the  polarization  equations 
p = tanh  x and  p = coth  x - l/x,  respectively.  The  broken  curve 
represents  the  Devonshire  approximation  to  the  J = 00  case. 

Fig . l8  - Comparison  of  the  spontaneous  polarization  of  triglycine  sulfate 
as  measured  by  Triebwasser  (Ref.  12)  with  the  prediction  of  the 
GMFT  using  the  Pade  approximant  form  of  the  local  field  and 


= 0.85 


p=  P/N/z 


Fig.  1 
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501  Commanding  Officer 

US  Army  Arctic  Test  Center 

002  APO,  Seattle  98733 

502  CO,  USA  Tropic  Test  Center 

ATTN:  STETC-MO-A  (Tech  Lib) 

Drawer  9h2 

001  Fort  Clayton,  Canal  Zone  09827 

503  Director 

US  Army  Adv-Matl  Concepts  Agcy 
ATTN:  AMXAM 

001  Washington,  . C.  20315 


501:  Commanding  General 

US  Amy  Materiel  Command 
ATTN:  AMCRD-R 

5001  Eisenhower  Ave. 

001  Alexandria,  VA  22333 


512  Commanding  Officer 

USACDC  Armor  Agency,  ATTN:  ATSAR-CTD 

001  Fort  Knox,  Kentucky  1:0121 

516  Commandant 

US  Army  Field  Artillery  School 
ATTN:  Target  Acquisition  Deot 

002  Fort  Rill,  Oklahoma  -73503 

517  Commanding  General 

US  Army  Missile  Command 
ATTN:  AMSMI-RFG  (Mr  N.  Bell) 

001  Redstone  Arsenal,  Ala.  35809 

518  Commanding  Officer 

Harry  Diamond  Laboratories 
ATTN:  AMXDO-RCB  (Pir.  Nemarich) 

001  Washington,  D.  C.  20li38 

596  Commanding  Officer 

US  Army  Southeastern  S’gnal  School 
A1TN:  ATSN-CTD-MS 
001  Fort  Gordon,  GA  30905 

598  Commanding  Officer 

USA  Satellite  Comm  Agency 
ATTN:  AMCPM-SC-3 
001  Fort  Monmouth,  N.  J.  07703 

599  Tri-Tac  Office 

ATTN:  CSS  (Dr.  Pritchard) 

001  Fort  Monmouth,  N.  J.  07703 
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607  Commanding  General 

USA  Tank -Automotive  Command 
ATTN:  AMSTA-Z,  Dr.  J.  Parks 

001  Warren,  Michigan  I4809O 


703  NASA  Sci  & Tech  Info  Facility 
ATTN:  Acquisitions  Br(S-AK/DL) 

P.0.  Box  33 

002  College  Park,  Maryland  207l»0 


610  Director 

Night  Vision  Lab  (USAECOM) 
ATTN:  AMSEL-NV- D 

001  Fort  Belvoir,  Virginia  22060 

6 lli  Chie  f 

Ofc  of  Missile  Electronic  Warfa 
Electronic  Warfare  Lab,  ECOM 
001  White  Sands  Missile  Range,  NM 
88002 

616  CG,  USA  Electronics  Command 
ATTN:  AMSEL-FP/P-IED 
(Mr.  C.  Mogavero) 

225  South  18th  Street 

001  Philadelphia,  PA.  19103 

617  Chief,  Intel  Materiel  Dev  & 
Ofc,  Electronic  Warfare  Lab 

001  Fort  Meade,  MD  20755 

680  Commanding  General 

US  Army  Electronics  Command 
000  Fort  Monmouth,  N . J . 07703 


705  Advisory  Gp  on  Electron  Devices 
201  Varick  St.  9th  Floor 

002  New  York,  New  York  1001 h 

706  Advisory  Gp  on  Electron  Devices 

ATTN:  Secy,  Sp  Gr  on  Opt  Masers 

201  Varick  Street 

002  New  York,  Nerw  York  lOOllj 

708  Ballistic  Msl  Radiation  Anal  Cen 
Univ  of  Mich.,  Willow  Run  Lab 
Institute  of  Science  & Tech 
001  PO  Box  6l8,  Ann  Arbor,  Michigan 
I|8107 

711  Mstals  and  Ceramics  Inf  Center 
Battelle 

Support  505  King  Avenue 
(ECOfQOl  Columbus,  Ohio  li3201 

712  Elec  Properties  Info  Center 
Hughes  Aircraft  Company 
Centinela  and  Teale  Streets 

001  Culver  City,  California  90230 


1AMSEL-NV-D 

1AMSEL-RD 

( eooM 

1AMSEL-NL-D 

J Ih-Iimm 

1AMSFL-WL-D 

1AMSEL-VL-D 

jjAMSEL-CT-D 

1AMSEL-BL-D 

1AMSEL-TL-D 

\ 

1AMSEL-TL-DT 

6 AMSEL-TL-ES 



(Ofc-  of  Record) 

1AMSEL-TE 

2AMSEL-KS-TI 
1AMSEL-GG-TD 
1AMSEL-EN 
/ :AMSEL-PA 

1USMC-LN0 

1 AMSEL-PP-C- 

ES-l 

715  Plastics  Tech  Eval  Center 

Picatinny  Arsenal,  Bldg  3ij01 

001  Dover,  N.  J.  07801  . 

717  Reliability  Analysis  Center 

Rome  Air  Development  Center 
ATTN:  J.M.  Schramp/RCRM 

001  Griffiss  AFB,  New  York  13MjO 

718  Remote  Area  Confligt  Info  Ctr 
Battelle  Memorial  Institute 
505  King  Avenue 

001  Columbus,  Ohio  ii3201 

719  Shock  and  Vibration  Info  Center 
Naval  Research  Lab  (Code  6020) 

001  Washington,  D.  C.  20390 


u 


* O*  nu**-*  euntrvci  Ada  COTR't 

mail  *ymbo* 


720  Thermophysical  Properties  Res  Ctr 
Purdue  Unlv.,  Research  Park 

2595  Yeager  Road 
001  Lafayette,  Indiana  47906 

721  Vela  Seismic  Info  Center 
University  of  Michigan 
Box  6 18 

003  Ann  Arbor,  Michigan  48107 


Institute  of  Science  6 Tech 
The  University  of  Michigan 
P.0.  Box  618  ( IR1A  Library) 

001  Ann  Arbor,  Ml  48107 

Director,  Night  Vision  Laboratory 
US  Army  Electronics  Command 
ATTN:  AMSEL-NV-FIR  (Dr.  R.  Ennulat) 
001  Fort  Bel  voir,  VA  22060 


Department  of  Electrical  Engineering 
ATTN:  Dr.  A.  van  der  Z 1 e 1 
University  of  Minnesota 

001  Minneapolis,  MN  55455  q 

Dr.  George  Hellmeler 
Director,  Advanced  Research 
Project  Agency 
1400  Wilson  Blvd. 

002  Arlington,  VA  22209 

C Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCAIC 

001  Fort  Leavenworth,  KS  66027 


Director,  Night  Vision  Laboratory 
US  Army  Electronics  Command 
ATTN:  AMSEL-NV-D  (Mr.  Edward  C.  Walker) 
Fort  Belvoir,  VA  22060 


Commander 

US  Army  Training  6-  Doctrine  Command 
ATTN:  ATCD-C1 
001  Fort  Monroe,  VA  23651 


D i rector 

Naval  Research  Laboratory 
CODE  2627 

001  Washington,  DC  20390 

Commander 
US  Army  R&0  Group 
(Far  East) 

003  AP0,  San  Francisco,  CA  96343 


Commandant 

US  Army  Engineer  School 
ATTN:  ATSEN-CTD 
001  Fort  Belvoir,  VA  22060 

Commanding  General 
US  Army  Materiel  Command 
ATTN:  AMCRD-TP  (Dr.  B.  Zarwyn) 
5001  Eisenhower  Ave. 

00!  Alexandria,  VA  22333 


5 


